Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 22 November 2012 



(MN WI^ style file v2.2) 



The nature of dark matter from the global high redshift HI 

21 cm signal 

M. Valdes^'^*, C. Evoli2'\ A. Mesinger\ A. Ferrara^'^ & N. Yoshida^'^ 

^ Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy 

^ II. Institut fiir Theoretische Physik, Universitdt Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany 
^ Kavli IPMU, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8568, Japan 

Department of Physics, University of Tokyo, 3-7-3 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 

22 November 2012 

ABSTRACT 

Wc study the imprint of dark matter (DM) annihilation on the global 21 cm signal 
from the Dark Ages to Cosmic Reionization. Motivated by recent observations, we 
■ focus on three DM candidates: (i) a 10 GeV Bino-like neutralino (ii) a 200 GeV Wino 

and (iii) a 1 TeV heavier particle annihilating into leptons. For each DM candidate 
we assume two values for the thermally averaged annihilation cross section (av), the 
standard thermal value {av)th = 3 x 10~^^cm'^s~^ and the maximum value allowed 
by WMAP7 data, {<Tv)max- We include the enhancement of DM annihilations due 
to collapsed structures, detailed estimates of energy deposition into the intergalactic 
medium (IGM), as well realistic prescriptions for astrophysical sources of UV and 
X-ray radiation. In these models, the additional heat input from DM annihilation 
suppresses the mean 21cm brightness temperature offset by STt ~ few-100 mK. In 
particular, the very deep STi, ^ —150 mK absorption feature at 20 ;$ z ;S 25 pre- 
dicted by popular models of the first galaxies is considerably reduced or totally erased 
by some of the considered DM candidates. Such an enhancement in IGM heating could 
come from either DM annihilations or a strongcr-than-cxpected astrophysical compo- 
nent (i.e. abundant early X-ray sources). However, wc find that the two signatures are 
not degenerate, since the DM heating is dominated by halos several orders of mag- 
nitude smaller than those hosting galaxies, whose fractional abundance evolves more 
slowly resulting in a smaller gradient: dSTh/dv ^ 4 mK/MHz in the range v^GO — SO 
MHz. The detection of such signals by future radio telescopes would be clear evidence 
of DM energy injection at high-redshifts. 

Key words: intergalactic medium - cosmology: theory - diffuse radiation - dark 
matter 



1 INTRODUCTION 

In the framework of the successful ACDM cosmology theory 
only ~ 5 % of the current energy content of the Universe 
is made of visible matter, while the rest is divided into two 
dark components, 73 % in the form of the so called Dark En- 
ergy, with the remaining 22 % accounted for by Dark Mat- 
ter (DM) which only interacts with baryons through grav- 
ity (Komatsu ct al. 2009). The existence of DM seems to 
be indirectly confirmed by a large set of observations (even 
though it is conceivable that in the future a new theory of 
gravity could do without DM), however its nature remains 
unknown. The accuracy of ACDM in explaining the evolu- 
tion of the Universe is one of the major successes of modern 
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cosmology; however the question of what makes up over 95 % 
of the Universe is still open. 

A huge effort has been made in the past decades to 
detect DM directly and indirectly. Currently, the best we 
can do is to place constraints on the properties (e.g. mass, 
annihilating cross section, decay rate) of some of the pro- 
posed DM candidates, thanks mostly to recent observations. 
We will focus here on Weakly Interacting Massive Particles 
(WIMPs) that are among the most popular DM candidates 
since their existence and properties are predicted by several 
extensions of the Standard Model of particle physics. 

WIMPs have a small but non negligible interaction 
cross-section with ordinary matter and therefore in principle 
they can be directly detected via elastic collision with the 
nuclei of terrestrial targets. Several recent experiments such 
as DAMA/LIBRA, DAMA/Nal, CDMS-II, EDELWEISS-II, 
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CoGeNT, XENONIOO have attempted such detection (see 
e.g. Bernabei et al. 2004, 2008, 2010; Aalseth et ah 2011a,b; 
Aprile et al. 2011; CDMS II CoUaboration et al. 2010; Ar- 
mengaud et al. 2010). Promising results came from Co- 
GeNT, which reported about a hundred events exceeding 
the expected background, possibly originated from the nu- 
clear recoil by scattering from DM particles of mass ^ 10 
GeV. A light WIMP with similar mass is also favored 
by the recent DAMA/LIBRA annual modulation signal 
(Bernabei et al. 2010; Aalseth et al. 2011a,b). 

Astrophysical observations in the X-ray and Gamma 
radiation bands have also been used to indirectly de- 
tect a DM generated signal. A lot of excitement followed 
the recent detection by the PAMELA, ATIC, FERMI- 
LAT and HESS experiments of an excess in the elec- 
tron/positron cosmic ray energy spectrum, which could 
be explained by annihilating DM with mass of order 
of 1 TeV (see e.g. Cirelli, Franceschini & Strumia 2008; 
Liu et al. 2010; Berg et al. 2009; Hooper & Tail 2009; 
Chen, Takahashi & Yanagida 2009; Abdo et al. 2010). The 
results appear to be far from conclusive due to the difficulties 
in properly modeling the cosmic ray energy spectrum and 
because other astrophysical sources such as pulsars could be 
responsible for the signal (see e.g. di Bernardo et al. 2011; 
Profumo 2012). 

Another promising observational window for the in- 
direct search of DM will be available in the near future 
when next generation radio interferometers such as the 
Low frequency Array (LOFAR^), the 21 Centimeter Array 
(21CMA2), the Murchison Widefield Array (MWA^) and the 
Square Kilometer Array (SKA*). These large interferome- 
ters will be able to map at arcminute scales the redshifted 
21 cm line corresponding to the hyperfine triplet-singlet line 
transition of the ground level neutral hydrogen (HI) at z > 6, 
i.e. during the Epoch of Reionization (EoR) and possibly 
well into the so called Cosmic Dark Ages of the Universe (see 
e.g. Peterson, Pen & Wu 2005; Bowman, Morales & Hewitt 
2006; Kassim et al. 2004; Wyithe, Loeb & Barnes 2005). 

The global sky-averaged 21 cm signal could poten- 
tially be measured as a function of frequency even by sin- 
gle dipoles: this is the main scientific aim of current in- 
struments such as the Experiment to Detect the Reion- 
ization Step (EDGES), while more advanced ones such as 
the Large-aperture Experiment to Detect the Dark Ages 
(LEDA^) and the Dark Ages Radio Experiment (DARE''), a 
lunar orbiting dipole experiment, are planned for construc- 
tion (see e.g. Bowman & Rogers 2010; Greenhill & Bernardi 
2012; Burns et al. 2012). 

Any release of energy from DM - for instance, as 
is the case considered in this work, through annihila- 
tions - would be in part absorbed by the intergalactic 
medium (IGM) as heat and ionization. If this oc- 
curred before the first astrophysical sources (z ^ 30), 
it would produce a deviation from the theoretically 
well established thermal and ionization history (see e.g. 

^ http://www.lofar.org/ 
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Chen & Kamionkowski 2004; Mapelli, Ferrara & Pierpaoli 
2006; Cirelli, locco & Panci 2009) and could alter signifi- 
cantly the HI 21 cm signal (see e.g. Shchekinov & Vasiliev 
2007; Furlanctto, Oh & Pierpaoli 2006; Valdes et al. 
2007; Natarajan & Schwarz 2009; Galli et al. 2011; 
Finkbeiner et al. 2012). If the DM energy injection 
occurred after the first astrophysical sources {z ^ 30), the 
deviation could be stronger and easier to observe, but 
would have to be disentangled from astrophysical uncer- 
tainties concerning the first galaxies. There are enormous 
observational challenges in detecting the 21 cm signal at 
z ^ 30, and such observations may be decades away from 
our current technological level. For this reason, we focus on 
the DM signal at lower redshift, after stars and galaxies 
started forming, and look for strategies to separate the DM 
signal from that of astrophysical sources. 

We investigate the effects produced on the HI 21 cm 
signal by three DM candidates: (i) a 10 GeV Bino-like 
neutralino (ii) a 200 GeV Wino and (iii) a 1 TeV heav- 
ier particle annihilating into leptons. These candidates have 
been recently proposed to explain the aforementioned indi- 
rect/direct hints of DM detection. We allow the annihila- 
tion cross-section to be in a range which is compatible with 
CMB observations. For the first time we compute the 21 
cm DM signal at z 30 in a self consistent scenario that 
includes a realistic prescription for the formation of astro- 
physical sources of UV and X-ray radiation. 

The details of the energy deposition into the IGM from 
DM annihilations are an essential ingredient to compute the 
HI 21 cm signal. In Evoli et al. (2012b) we have studied 
the cascade produced by the products of DM annihilations 
for the three considered DM candidates: this allows us to 
discuss the impact of DM annihilations in a physically con- 
sistent way and to assess the observability of the chosen DM 
candidates with HI 21 cm observations in the near future. 

The rest of the paper is organized as follows. In Sec. 2 we 
briefly introduce the DM candidates which we study; in Sec. 
3 we explain in detail the method that we follow to compute 
the 21 cm DM signal. In Sec. 4 we present the results of our 
calculations, and in Sec. 5 we discuss the results and draw 
our conclusions. 



2 DARK MATTER MODELS 

In this Section we briefly introduce our DM models of choice. 
For a more detailed discussion of such DM particle candi- 
dates we refer the reader to Evoli et al. (2012a), were they 
are described extensively. We select three sample cases which 
have been recently investigated in connection with hints of 
DM signals in either direct and indirect search experiments. 
These sample cases are also representative of three different 
WIMP mass regimes, ranging from fairly light models to 
multi-TeV DM, and of three different kinds of annihilation 
channels. 

Wino: We consider a pure Wino within the Minimal 
Supersymmetric extension to the Standard Model (MSSM). 
The recent interest in this model has been stimulated, be- 
sides its peculiar signatures at the LHC (Bertone et al. 
2012), by the claim (Grajek et al. 2009; Kane, Lu & Watson 
2009) that a Wino with mass of about 200 GeV can explain 
the rise detected by PAMELA in the positron fraction. 
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Bino: We consider the bh case with mass of 10 GeV to 
model a strong couphng with quarks as suggested by recent 
results of direct detection experiments (e.g. CDMS, DAMA). 

Leptophilic: Again motivated by the PAMELA 
positron excess, and possibly in connection with the local all- 
electron (namely electrons plus positrons) flux measured by 
Fermi and HESS, several analyses have considered the possi- 
bility of very heavy dark matter WIMPs, with masses up to 
several TeV and very large pair annihilation cross section 
(see e.g. Cirelli et al. 2009; Bergstrom, Edsjo & Zaharijas 
2009). The results of such studies are that, to account for 
the electron/positron component without violating the an- 
tiproton bounds, dark matter needs to be leptophilic, i.e. 
the final products of the annihilation being dominantly lep- 
tons, most likely a combination of e^e~ and jj.'^ fi~ . It has 
been recently pointed out (Ciafaloni ct al. 2011) that for 
very heavy WIMPs the radiative emission of soft electroweak 
gauge bosons is crucial to model the high-energy spectrum. 



3 METHOD 

We first constrain the properties of our DM candidates from 
CMB data. This allows us to associate to each particle a 
range of allowed values for the annihilation cross section 
(cTv). We then use the fractional energy depositions from 
MEDEA2 (Evoli et al. 2012b) and the galactic radiation fields 
from 21cmFAST (Mesinger, Furlanetto & Cen 2011) in order 
to compute the thermal and ionization evolution of the IGM 
and the associated global 21cm signal. We describe each of 
these steps in turn. 

3.1 DM heating 

Including the formation of substructures at redshift 50 
naturally enhances the effects of DM annihilations: while 
they do not dramatically impact the global reionization 
history (Cirelli, locco & Panel 2009) the increased energy 
injection could alter predictions of the observable 21 cm 
signal by heating the IGM. The total energy release by 
DM annihilations per unit volume is given by (see e.g. 
CireUi, locco & Panel 2009; Chluba 2010): 

^ = 2M^c\cjv)nlil + B{z)), (1) 

where (av) is the thermally averaged annihilation cross sec- 
tion, is the DM particle mass and = no ,x(l + is 
the number density of DM particles and anti-particles with 
present day value: 

no,. = 1.2xlO «cm . (2) 

The term B{z) defines the effective averaged DM density 
resulting from structure formation and can be written as: 

B{z) = ^^r dMM^{z,M)F{M,z) (3) 

where An/dM is the halo mass function obtained adopt- 
ing the Press-Schechter formalism (Press & Schechter 1974) 
and F{M, z) is the concentration function depending on the 
distribution of DM inside halos, a function that peaks typ- 
ically for sub-halo masses Msh ~ 1 M0 at z = 20 (see 



Fig. 7; we refer the reader to Cirelli, locco & Panel 2009 
for details)^. We use here a prescription for the effect of 
substructures and consider a minimum mass for primordial 
proto-haloes of Afmin = 10~® M0 for the 10 GeV Bino, and 
A/min = 10"^ M0 for the other two more massive candi- 
dates. This choice is motivated by recent calculations of the 
exponential cutoff mass in the power spectrum as a function 
of neutralino mass, resulting from free-streaming of the DM 
particles after the kinetic decoupling (see e.g.Bringmann 
2009; van den Aarssen, Bringmann & Goedecke 2012). 

Recent works (e.g. Slatyer, Padmanabhan & Finkbeiner 
2009; Galli et al. 2009) have investigated the possible role 
of Sommerfeld enhancements due to which (ctu) can be 
considerably boosted and becomes a function of redshift. 
We choose not to make assumptions on the physical 
processes that could boost {av) and treat it as a parameter 
constrained by CMB observations. 

We covered until now the energy production by DM an- 
nihilations. The most relevant quantity for our purposes is 
however the energy deposition into the IGM. Only a small 
fraction of the energy released by DM annihilations is finally 
deposited into the IGM in the form of heat, excitations and 
ionizations either of hydrogen and helium. The absorbed 
fraction depends on the DM candidate - in the form of the 
initial energy spectrum from the annihilation event - and on 
the environment where the annihilation takes place, specifi- 
cally on the ionized fraction of the ambient gas and on the 
energy density of CMB photons, important for the Inverse 
Compton (IC) scatterings. 

To model this, we use the Monte Carlo scheme MEDEA2 
(for details see Evoli et al. 2012b). Through MEDEA2 we 
are able to calculate what fraction of the energy released 
by a single annihilation event goes into: (i) H and He 
ionizations (ii) excitations (ill) heating. A set of handy 
fitting formulae that take into account the dependence 
on z and 2:0 are then given for the respective quantities 

fio?i{Xe,z), fa{Xe,z), fh{Xe,z). 

3.2 CMB constraints 

We compute the modifications induced by DM annihilations 
on the CMB power spectra and on the integrated Thomson 
optical depth to verify that our models are consistent with 
WMAP7 observations and to put upper limits on (av) for 
the considered DM candidates. To do so we modify the pub- 
lic code CAMB (Lewis & ChaUinor 2011) and add a DM term 
to the evolution equations of the IGM kinetic temperature 
and ionized fraction (Fig.l). 

With this modified version of CAMB code we can cal- 
culate the effects of DM annihilations on the temperature 
(TT), polarization (EE) and temperature-polarization (TE) 
CMB spectra. The TT spectrum allows us to put the most 
sensitive constraints of DM properties (see Fig.2). Our re- 
sults are obtained assuming that the cosmological parame- 
ters have best-fitting values as indicated by the 7-yr WMAP 

^ Although extrapolating the halo mass function down to 
such small masses is highly uncertain, wc note that at the 
scales and rcdshifts of interest, the Press-Schechter and Sheth- 
Thormen cumulative mass functions (Press & Schechter 1974; 
Sheth & Tormen 1999; Jenkins et al. 2001) agree to ~ 10%. 
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Table 1. The DM models described in Sec. 2. 
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Figure 1. IGM ionization fraction as a function of rcdshift. Black 
solid line is the case without annihilating DM, for the color scheme 
of the DM models see Tab. 1. 



data (Komatsu et al. 2011). For each DM candidate we in- 
crease the value of the thermally averaged annihilation cross 
section from the standard value {av)th = 3 x 10~'^®cm^s~^ 
up to a value (av)m.ax for which the computed TT CMB 
spectrum exceeds by more than 3o" from the best-fit the 
WMAP7 data. The upper limits {ijv)max for each of the con- 
sidered DM candidates are given in Tab. 1 along with the 
color and line-style later used in the plots. CMB observa- 
tions by Planck will allow to put more stringent constraints 
on the properties of the DM candidates in the near future 
(The Planck Collaboration 2006; Galh et al. 2009). 

We then calculate the contribution of DM to the Thom- 
son optical depth, 5re, by integration: 



— dz 

az 



(4) 



where 5xe. is the difference in the ionized fraction between 
the standard scenario and the case in which DM annihi- 
lates. We assume here that the Universe is fully ionized at 
redshift z < 6: this contributes by a factor Te ^ 0.04 to 
the WMAP measured total optical depth = 0.088 ±0.015 
(Komatsu et al. 2011; ShuU et al. 2012). Note that DM an- 
nihilations are only imprinted on the CMB as an additional 
source of ionization. Therefore the above procedure is ef- 
fectively the same as choosing an upper limit for the extra 
contribution to Te from DM annihilations. From Tab. 1 we 
see that this is Ste ~ 0.06, making the total optical depth 
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Figure 2. (Top panel) TT CMB power-spectrum for the DM 
models in Tab. 1. Black solid line is the case without annihilat- 
ing DM. (Bottom) Residual of the CMB power-spectrum with 
respect to the case without DM annihilations. The DM models 
are presented in red (10 GeV); blue (200 GeV); green (1 TeV). 
The triangles and squares are used for {av)tfi and (av) 

max re- 
spectively. 



(in the absence of any additional astrophysical sources at 
2 > 6): Te ~ 0.1. This limit is conservatively low, roughly 
corresponding to the 1-a upper limit from WMAP7 obtained 
from the TE cross-correlation (Komatsu ct al. 2011). 

Notice that the heating/ionization contribution of the 
two heavier DM candidates assuming (av) max increases pro- 
portionally by a larger factor over the case with {av)th than 
what found for the 10 GeV Bino when assuming {(Tv),nax 
rather than {av)th- This is due to our assumption of a 
smaller minimum halo mass Mmin ~ 10"^ Mq that enhances 
strongly the effects of substructures and to which the CMB 
constraints are, on the other hand, less sensitive. 

3.3 21 cm signal 

One of the observable quantities most likely to carry a trace 
of the effects of DM annihilations is the redshifted 21 cm line 
associated with the hyperfine transition between the triplet 
and the singlet levels of the neutral hydrogen ground state. 
This signal is most commonly expressed in terms of the dif- 
ferential brightness temperature between a neutral hydrogen 
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patch and the CMB (neglecting redshift-space distortions) 
Ts — TcMB , 



5n 



l + z 



-(1 



272;hi(1 + ^) 1- 



^c^ 



Ts 



l + z 0.15 
10 ^Mh? 



1/2 



0.023 



mK, 



(5) 



where xm is the neutral fraction of the gas, 5 the overdensity, 
and Ts is the spin temperature which is set by the number 
densities of hydrogen atoms in the singlet (no) and triplet 
(ni) ground hyperfine levels, ni/no = 3 exp(— O.OGSK/Ts). 

It is theoretically well known that in the presence of the 
CMB alone, Ts reaches thermal equilibrium with Tcmb = 
2.73 (1 + z) K on a short time-scale, making the HI unde- 
tectable in emission or absorption. However, collisions and 
scattering of Lya photons — the so-called Wouthuysen-Field 
process or Lya pumping (e.g. Wouthuyscn 1952; Field 1959; 
Hirata 2006) can couple Ts to Tk- 

The spin temperature can be written as (e.g. 
Furlanetto, Oh & Briggs 2006): 



m -1 Tcmb ^ + XqTq ^ + XcTk ^ 
-'S ~ — ; 



1 ~\~ Xq, ~\' Xc 



(6) 



where Ta is the color temperature, which is closely coupled 
to Tk (Field 1959), and Xa and Xc are the two coupling 
coefficients relative to Lya scattering and collisions respec- 
tively. If either collisions or Lya radiation couple Ts to Tk 
then the neutral hydrogen will be visible in absorption or 
emission depending on whether the gas is colder or hotter 
that the CMB respectively. 

For details about the physics associated with the 
HI 21 cm line and with the determination of Ts, 
the quantity that governs it, we refer the reader 
to, e.g., Furlanetto, Oh & Briggs (2006); Hirata (2006); 
Valdes & Ferrara (2008). For our purposes it is important 
to state here that the physical quantities that determine Ts 
are: (i) the gas density uh', (ii) the CMB temperature Tcmb; 
(iii) the kinetic temperature of the gas, Tk; (iv) the ionized 
fraction Xe; and (v) the Lya background intensity Jq. 

While we know the average evolution of uh and Tcmb 
we have to determine the others as a function of red- 
shift. The equations that describe the average evolution 
of the ionized fraction Xe and of the kinetic temperature 
Tk are the following (see e.g. Chen & Kamionkowski 2004; 
Mesinger, Furlanetto & Cen 2011): 



dTK 
dz 



dXe 

dz 
2Tk 



dt 
dz 



Fio 



- aBSeJlfc/Hl 



(7) 



+ 



dt 



l + z Sksil + fHe + Xe) dz 



where rib ~ nb,oil + z)^ is the mean baryon number density, 
ep{z) is the heating rate per baryon for each process p in erg 
s~^, Fion is the ionization rate per baryon, ob is the case-B 
recombination coefficient, fcs is the Boltzmann constant and 
/hb is the helium fraction by mass. 

The term Fion includes both the contribution from 
galaxies and the term that accounts for DM annihilations. 
Similarly we have that ep{z) = ecMB(z) + ex(z) + eDM(z) is 



the sum of three contributions: (i) ecMB (z) is Compton heat- 
ing from CMB photons; (ii) ex(z) is heating from astrophys- 
ical sources, which we take to be dominated by X-rays; (iii) 
eDM(z) is DM heating. Notice that when including DM in 
the terms Fion and ep{z) we use the specific fractional energy 
depositions fh{xe,z), fion{xe,z) from Evoli et al. (2012b). 

The last equation needed to compute the 21 cm back- 
ground is the one describing the evolution of the Lya back- 
ground intensity Jc,: 



Ta — Ja,R ~\~ 'Ja,C ~\~ Ja 



Ja.X + Ja 



(8) 



where the contributions on the RHS correspond to recom- 
binations, coUisional excitations by electron impacts, direct 
stellar emission, X-ray excitation of HI, and DM annihila- 
tions (respectively) (see e.g. Madau, Meiksin & Rees 1997; 
Valdes et al. 2007; Mesinger, Furlanetto & Cen 2011). The 
coupling coefficient, Xa is proportional to the Lyman a back- 
ground flux, Ja (e.g. Furlanetto, Oh & Briggs 2006). 

In general, if we neglect the energy input from DM, we 
expect Tk and Ts to track Tcmb down to 2 ~ 300, when 
Tk decouples from Tcmb and starts decreasing adiabatically 
as the Universe expands. Ts is then coupled to Tk due to 
the high gas density and the consequent strong coUisional 
coupling. At z ~ 70 Ts gradually couples to Tcmb until, at 
2 ~ 30 becoming virtually identical to it. It is believed how- 
ever that at around this redshift the first collapsed luminous 
sources would ignite. Radiation produced by the first galax- 
ies starts to ionize and heat the gas: on a short timescale, 
at a redshift zwf 25 — 30, the Lya pumping effectively 
couples Ts to Tk, and only later, at Zhoat ~ 18 — 22, heating 
from galaxies drives Tk to values higher than Tcmb , making 
the neutral regions in the IGM visible in emission. Therefore 
a second absorption feature is expected at Zhoat ^ z ^ zwf 
(see e.g. Furlanetto, Oh & Briggs 2006; Pritchard & Loeb 
2008). We can divide the global evolution of 5Tb in six main 
phases: 



(i) STb = for z > 300; 

(ii) 5Tb < for 30 ^ z 300: this extended absorption 
feature has a minimum of ~ 45 mK at z ~ 90; 

(iii) 5Tb ~ for zwp ;S z < 30; 

(iv) 5Tb < at Zhoat ^ z ^ Zwf, a second absorption fea- 
ture less extended in redshift than the previous one but 
deeper, with 5Tb,,nin ~ —150 mK; 

(v) 5Tb > 0, i.e. in omission, due to heating by galaxies at 



down to the Epoch of Reionization at z — zeor ~ 

6-8; 

(vi) 5Tb ~ for ^ z zeoR, where we still have a signal 
only from self-shielded systems. 

Introducing the effects of DM annihilations in the stan- 
dard scenario described above can produce deviations on 
5Tb: in particular the second absorption feature can be 
strongly modified by energy release by DM annihilations. 



at a redshift range 15 



25 that will be probed by the 



next generation of radio observatories such as LEDA and 
SKA. 
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Figure 3. Slices through the 21cm brightness temperature maps for two models with different prescriptions for the galaxy properties 
(Mesinger, Ferrara & Spiegel 2012). The simulations are 750 Mpc on a side, with a resolution of 500"^ . Each slice is 1 cell (1.5 Mpc) 
thick. The horizontal axis shows evolution along the comoving linc-of-sight coordinate. The top panel corresponds to a "fiducial" model, 
in which the X-ray luminosity of primordial galaxies is the same as that observed in nearby starburst galaxies, while the lower panel 
corresponds to a model in which primordial galaxies arc much more efficient in generating X-rays, saturating the soft X-ray background 
at 2: 10 (see text for details). 



3.4 Radiation from astrophysical sources 

To compute the astrophysical contribution to Fion, fp, and 
Ja, we use the publicly available code, 21cmFAST This code 
uses a combination of perturbation theory and excursion-set 
formalism to compute various cosmic fields, and is in good 
agreement with radiative transfer simulations of reioniza- 
tion (e.g. Mesinger & Furlanetto 2007; Zahn et al. 2011; 
Mesinger, Furlanetto & Cen 2011). The code is fully de- 
scribed in Mesinger, Furlanetto & Con (2011), to which we 
refer the interested reader for more details. Here we briefly 
note that the code takes into account inhomogeneous X-ray 
ionization and heating, as well as Lya pumping from the 
first UV sources, integrating the evolution of cosmic struc- 
tures and radiation fields along the light cone. Although 
21cmFAST computes 3D realizations, in this work we only 
study the global 21cm signal, deferring analysis of the spa- 
tial structure to future work. 

The simulation boxes are 750 Mpc on a side, with a 
final resolution of 500"^. The initial conditions are sampled 
on a 1500^ grid. In Fig. 3 are presented the 1-cell thick (1.5 
Mpc deep) slices through the HI 21 cm brightness tempera- 
ture maps for two models with different prescriptions for 
the galaxy properties (Mesinger, Ferrara & Spiegel 2012). 
The top panel corresponds to a "fiducial"^ model, in which 
the X-ray luminosity of primordial galaxies is the same as 

® http : //homepage . sns . it/mesinger /Sim. html 
^ Notice that our "fiducial" model does not exactly correspond 
to the fiducial model in Mesinger, Ferrara & Spiegel 2012 (consis- 
tent with the recent measurement of the 0.5 — 8 keV X-ray lumi- 
nosity of star forming galaxies from Minoo, Gilfanov & Sunyaev 
2012), but instead corresponds to their Tle4_fuvl_fx5_lkeV 
model (based on an extrapolation of the 2 — 10 keV data from 
Gilfanov, Grimm & Sunyaev 2004). This model represents a more 
conservative choice since the astrophysical heating is 4 times 
higher, thus making the DM annihilation signal less apparent. 



that observed in nearby starburst galaxies (e.g. Furlanetto 
2006 and references therein). The lower panel corresponds 
to an "extreme" model in which primordial galaxies, al- 
beit rarer and appearing later, were much more efficient in 
generating hard X-rays. The later scenario is inspired by 
recent theoretical (e.g. Linden, Profumo & Anderson 2010; 
Mirabel et al. 2011) and observational (Reichardt et al. 
2011; Mesinger, McQuinn & Spergel 2012) claims. More 
specifically, these two models have the following character- 
istics: 

(i) fiducial (top panel): Galaxies hosting UV and X-ray 
sources reside in atomically-cooled halos with virial tem- 
peratures Tvir > 10* K (corresponding to halo masses of 
A/halo > 3 X lO^M© at 2 20). Lya pumping (i.e. spin 
temperature coupling) is dominated by early UV sources, 
assuming PopII stellar spectra (e.g. Barkana & Loeb 
2005) and a 10% efficiency of conversion of gas into 
stars. The X-ray luminosity of galaxies follows a hv'^'^ 
power law shape with a lower limit of hi/o ~ 300 eV (e.g. 
Madau et al. 2004), and an X-ray efficiency corresponding 
to ~2 X-ray photons per stellar baryon. Similar mod- 
els, inspired by lower-redshift X-ray binary-dominated 
star-burst galaxies, have been explored in prior analytic 
studies (e.g. Furlanetto 2006; Pritchard & Furlanetto 2007). 

(ii) extreme (bottom panel): Galaxies hosting UV and 
X-ray sources reside in more massive halos, similar to 
ones observed at moderate redshifts (e.g. Labbc ct al. 
2010) and which are more resilient to feedback effects 
(e.g. Springel & Hernquist 2003; Mesinger & Dijkstra 2008; 
Okamoto, Gao & Theuns 2008). These galaxies have virial 
temperatures Tvir > 10^ K, corresponding to halo masses 
of Mhaio > 10^ Mq at z « 20. Interestingly, in this case 
Lya pumping is dominated by the X-ray excitation of HI, 
and has a very different spatial signature from the fidu- 
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cial model. The X-ray luminosity of galaxies also follows 
a hu^^'^ power law shape, but with a more energetic lower 
limit of hvo — 900 eV (corresponding to heavy obscura- 
tion, e.g. Lutovinov ct al. 2005), and an X-ray efficiency 
corresponding to ~4000 X-ray photons per stellar baryon. 



This scenario is considered "extreme" since the 



10 X- 



ray sources dominate reionization and saturate the unre- 
solved soft X-ray background (Hickox & Markevitch 2007; 
Mcsingcr, Fcrrara & Spiegel 2012). 

Both scenarios are consistent with the WMAP7 constraints 
on Te at 2a (Komatsu et al. 2011). Notice that the values of 
{crv)max computed in Section 3.2 represent an upper limit 
and are not self-consistent with the contribution to Te from 
the astrophysical sources modelled by 21cmFAST . On the 
other hand the DM models with the thermal cross-section 
{(jv)th are fully consistent since they only contribute negli- 
gibly to Te. 

It is obvious from Fig. 3 that the fundamental epochs 
in cosmic evolution show very different 21cm signatures in 
the two models: spin temperature (WF) coupling (black ^■ 
yellow); X-ray heating of the IGM (yellow — ^ blue); reioniza- 
tion (blue — >■ black). Here we focus on the global HI 21 cm 
signal as a way to constrain DM and save spatial signature 
for future work. 



4 RESULTS 

We present here the results obtained for the three DM mod- 
els introduced in Sec. 2, each with two values for the an- 
nihilation cross (av): the standard thermal value, and the 
maximum allowed by CMB constraints as described in Sec 
3.2. We solve the evolution equations of the kinetic tem- 
perature and ionized fraction of the IGM for each of these 
DM candidates using the fractional energy depositions from 
the code MEDEA2 and use a prescription for the formation 
of collapsed sources of radiation from the publicly available 
code 21cmFAST , that allows us to include consistently X-ray 
heating, ionization and Lya pumping from galaxies. 

In Fig. 4 we compare the effects produced on Ts by 
our annihilating DM candidates. To keep the figure simple 
we don't plot the modified Tk curves, corresponding to the 
annihilating DM cases. However the behavior is simple: the 
kinetic temperature at high redshift increases because of the 
high density of DM and the corresponding higher chance 
of annihilations. As the Universe expands and its content 
is diluted Tk simply settles on a slightly higher adiabat. 
When the Lya pumping becomes efficient a.t z ^ zwf, then 
Ts perfectly tracks Tk, which increases sharply at z ^ -Zhoat- 

The additional heating from annihilating DM particles 
drives Ts closer to Tcmb for 40 < z < 100. Therefore we can 
expect a reduction of the absorption feature of the homo- 
geneous 21 cm background in that redshift range. A larger 
thermally averaged cross-section corresponds to stronger ef- 
fects on Ts for an annihilating DM candidate of a given 
mass. At the same time, heavier DM candidates deposit their 
energy into the IGM less efficiently (see Evoli ct al. 2012b), 
therefore given a certain value of (av) the lower mass DM 
candidate will generally produce larger deviations on Tg. 
The effects are much more evident at z ^ 30: galaxies start 
to form and quickly drive Ja to values high enough for effi- 
cient Lya pumping. As a consequence Ts decreases sharply 




Figure 4. Tg, Tj, (thin black solid line) and Tcmb (black short- 
dashed line) as a function of redshift. The colored Tg curves show 
the modified behavior due to DM annihilations following the color 
scheme given in Tab. 1. The solid thick black line is the standard 
Ts without DM energy input. The top and bottom panels repre- 
sent the calculations done for the fiducial and extreme models re- 
spectively. The box on the lower right corner in each panel shows 
the behavior of Ts for the considered DM candidates without 
including radiation from luminous sources. 



tracking Tk at a time when the combined heating from DM 
annihilations and from galaxies is still not sufficient to heat 
the gas above Tcmb- Eventually the same heating sources 
make Tk ~ Ts > Tcmb. The different models in the upper 
and lower panels of Fig. 4 show different behaviors since in 
the "extreme" case Lya coupling is achieved later, and the 
heating starts earlier, making the region at 16 ^: 2 30 in 
which Ts < Tcmb both shallower (in K) and narrower (in 
redshift). 

This is reflected directly on the behavior of STt, shown 
in the panels of Fig. 5, for which we give a quantitative 
analysis case by case. A summary of the results is available 
in Table 2. 



4.1 10 GeV Bino 

The lightest considered DM candidate, the 10 GeV Bino, 
produces the largest signature in the "flducial" model (top 
panel of Fig. 5), both when assuming a standard ther- 
mal cross section {av)th = 3 x lO^^^cm'^s^^ (solid red 
curves in Fig. 5) and when taking into account the max- 
imum cross section allowed by CMB data, in this case 



1.0 X 10 ^^cm^s ^ (dashed red curves). In the 
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first case the signal is 5Tt ~ 4 — 10 mK on a redshift range 
45 ^ z ^ 300, with a peak of ^ 10 mK ai z ^ 100, and is 
much more substantial for the second absorption feature at 
16 ^ 2 ^ 30, with a deviation with respect to the standard 
Ccise in which DM does not annihilate (which we denote here- 
after STbfi) of ATi,,DM = \6Tb - 5Tb,o\ ~ 100 mK, a large 
signal at frequencies ~ 80 MHz. The effect is enhanced 
essentially by a factor two for the higher allowed annihila- 
tion cross section and reaches values of ATb^uM ~ 20 mK at 
45 ^ z ^ 300, while the second absorption feature at lower 
redshift is essentially erased, with the IGM appearing in 
even emission already by z ^ 25. These very large signals, 
both for {av)max and for {av)th could be detected by future 
radio observations (see Section 4.5). In the "extreme" model 
case (bottom panel) the DM signature before the first astro- 
physical sources turn on {z ~ 45 — 300) is obviously identical 
to the "fiducial" model case. The second absorption feature 
instead changes substantially and is both shallower, with a 
minimum value of 5Tb ~ —60 mK at z --^ 21, and narrower 
in redshift at 2: ~ 18 — 25. This reflects on the DM signal: 
when assuming {av)th {{crv)max), ATt^DM ~ 30 (45) mK at 
z-21. 

4.2 200 GeV Wino 

The case of the 200 GeV Wino, described by the blue curves 
in Fig. 5 produces less evident deviations: assuming {av)th 
we can see no effect on the STb as the solid green line is 
virtually coincident with the solid black line for both reion- 
ization models. Obviously this translates into a constant 
ATb,DM ~ 0, with no chances of direct detection. On the 
other hand the case {av)max = 1-2 x 10~^*cm^s~^ shows a 
ATb,DM = 2 — 9 mK for 60 2: ^ 300 and a massive devia- 
tion at 2: ^ 23 that practically erases the second absorption 
feature and drives the signal to emission already at 2; ~ 20. 
This behavior is present for both the fiducial and extreme 
reionization models. 

4.3 1 TeV Leptophilic 

Our most massive candidate, the heavy DM particle of rest 
mass 1 TeV that pair annihilates into leptons, is the one 
that produces by far the smallest deviations on the HI 21 cm 
background: for {av)th deviations are negligible while for the 
maximum allowed cross section, in this case {av),nax = 1.4 x 
10~^^cm'^ s~^, they reach ATb^oM ~ 15 mK in the range 
60 <z < 300 and ATb,DM ~ 65 (25) mK at 20 < 2 < 30 for 
the fiducial (extreme) reionization models. 

As mentioned previously, the fractional increase in the 
DM signal for the two heavier candidates assuming {av)max 
rather than {av}th is larger than what found for the 10 GeV 
Bino, due to our choice for the heavier DM particles of a 
smaller minimum sub-halo mass A'/n.in = lO-'' Mq. 

4.4 Isolating the DM annihilation signal 

4-4- 1 Degeneracy with astrophysics 

Although the additional heat input from DM annihilation 
is easy to isolate in the well-understood epoch of the Dark 
Ages (before the first astrophysical sources, z Si 30-50), the 
Earth's ionosphere makes observations of such high redshifts 
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Figure 5. STb ^ a- function of redsliift for all the considered DM 
models. The standard case with no energy input from DM annihi- 
lations is represented with the thick black line. The colored lines 
follow the scheme given in Tab. 1. The top and bottom panels rep- 
resent the calculations done for the fiducial and extreme models 
respectively. The red dotted curve in the upper panel corresponds 
to the thermal Bino model, but ignoring the X-rays from astro- 
physical sources (i.e. assuming astrophysical sources of X-rays 
ignited at later times). Notice that some of the curves relative to 
the DM models are below the black curves at z = 30 — 40. This 
is due to the term Ja,DKU i-c the extra Lya coupling from the 
DM annihilations. 

very challenging (though there are a couple of non-terrestrial 
radio telescopes being considered, such as DARE and the 
Lunar Radio Array^", LRA). The lower-redshift signal at 
Zheat < z < zwF is easicr to observe; however in this case 
the DM annihilation heating must be disentangled from un- 
certainties in astrophysics. This is evident when compar- 
ing the top and bottom panels of Fig. 5. In the "extreme" 
model in the lower panel the emission from astrophysical 
sources produces a much shallower absorption feature, and 
the black solid line could be confused with the red solid 
line (corresponding to the case of 10 GeV Bino annihila- 
tions with {av)th) in the upper panel. Therefore detecting 
a ~ —50 mK global absorption signal at z ~ 22 could be 
an indication of 10 GeV DM annihilations or of a strong X- 
ray emission by primordial galaxies. Notice however that the 
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X-ray emission can be constrained in several ways, e.g. by 
tighter constraints on the unresolved soft X-ray background 
of or by a study of the topology of the HIT regions once the 
future radio interferometers start probing the tomography 
of the IGM beyond the EoR. Unfortunately, both of these 
observations are only indirect. The soft X-ray background 
at z = only constrains the high-energy component of the 
source's spectrum (~ 20-40 keV at z ~ 20), which does not 
interact with the IGM (e.g. McQuinn 2012). Likewise, the 
EoR at 2 ~ 10 could probe a different source population 
than was present at z ~ 20 (e.g. Ricotti & Ostrikcr 2004). 
Nonetheless, if an upper limit is given to the efficiency of 
the X-ray emission by primordial galaxies a strong reduc- 
tion of the depth of the absorption feature would be a clear 
indication of DM annihilations. 

4-4-2 DM signal gradient 

It is interesting to notice that purely by coincidence galax- 
ies and DM annihilations start to heat the IGM at about 
the same redshift in these models even though the pro- 
cesses are entirely different and driven by structures with 
different mass. As previously mentioned, galaxies with mass 
Mhaio ~ 3 X 10^ Mq (10^ Mq) are responsible for most 
of the astrophysical X-ray heating at z = 20 for the "fidu- 
cial" ("extreme") reionization model. On the other hand 
the function F{M, z) in Eq. 3 peaks strongly for substruc- 
tures Mah ~ 1 M© at the same redshift, as we show in Fig. 
7. The difference in the collapsed fraction above these dis- 
parate mass scales is made up by the coefficients in the heat- 
ing rates. In Fig. 6 we present the heating rates per baryon 
for 10 GeV DM annihilations (edm) and for galactic X-ray 
heating (ex) for the fiducial reionization model, compared 
with the adiabatic cooling rate (see Eq. 7). Notice that the 
heating rate from astrophysical sources of X-rays is much 
steeper than the DM heating, but they become dominant 
over the adiabatic cooling at about the same redshift. We 
also plot the heating rate relative to 200 GeV DM annihi- 
lations to show that the slope is similar for different DM 
candidates. 

These differences in the slope are easy to understand: 
the fractional increase of the collapsed fraction in ^ IM0 
halos which drive the DM heating is much slower than the 
fractional increase in the high-end tail of the mass function 
(i.e. the halos which host the first galaxies). This difference 
is fundamental, and presents an unambiguous way of disen- 
tangling heating from astrophysics and heating from DM. 

We quantify this further with the dotted red line in the 
upper panel of Fig. 5. This curve corresponds to the case 
of annihilating 10 GeV Binos with {av)th, keeping the Lya 
pumping due to stellar sources but switching-off astrophys- 
ical X-ray heating, which could indeed take place at lower 
redshifts. The increase of STt aX z £j22 {v 60 MHz) here 
is only due to heating from DM annihilations; the different 
slope could be a clear indication of DM heating. 

In the bottom panel of Fig. 8 we further investigate 
the difference in slope between the DM and galactic heat- 
ing and compare them with the standard scenario in which 
DM does not annihilate. We study the gradient of the signal 
with frequency, i.e. dSTb/dv in mK/MHz, for the 10 GeV 
Bino, assuming the standard thermal annihilation cross and 
the fiducial reionization model, and keeping (solid line) or 
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Figure 6. Cooling/heating rates as a function of redshift. We 
compare here: (i) heating rate et{z) from galaxies for the fiducial 
reionization model (black dotted line); (ii) eDM(z) for 10 GeV 
DM annihilations assuming a standard {crv)th with and without 
taking into account the boost from substructures (dot-dashed and 
dashed red lines respectively); (iii) heating from 200 GeV Wino 
annihilations including the effect of substructures and assuming 
{'^^')th (dot-dashed blue line); (iiii) the adiabatic cooling rate of 
the expanding gas (solid line). Compton heating is negligible at 
these redshifts. By coincidence heating from 10 GeV DM annihila- 
tions and from galaxies becomes dominant over adiabatic cooling 
at a similar redshift, z ~ 20. 



switching off (dotted line) the heating from astrophysical 
sources. We choose the 10 GeV particle because it is the 
most promising candidate since it produces a strong signal 
at 2 ~ 20 with {crv)th, and we assume the "fiducial" model 
because in the "extreme" case the X-rays drive the heating 
and ionization of the IGM and produce with secondary in- 
teractions the Lya flux that couples Ts to Tk, therefore it 
would be incorrect to keep the Lya pumping and switch off 
the heating since they are driven by the same sources. 

It is evident from the figure that the inclusion of DM 
heating confines the slope to dSTt/dv ^ 4 mK/MHz (red 
curves) , whereas without DM heating the slope reaches val- 
ues of dSTt/dv ~ 14 mK/MHz. This slope would only in- 
crease if rarer (more massive) galaxies drive X-ray heating 
(such as in the "extreme" model). The detection of such a 
small peak gradient (dSTb/dv ~ few mK/MHz) could provide 
clear evidence of DM annihilation. 

Including the effects of star formation in smaller sized 
mini-halos would not alter things significantly, since the frac- 
tional increase of the collapsed fraction in mini-halos is still 
much steeper than that relative to the Si IM© sub-halos 
responsible for the DM heating. In principle, a scenario 
in which the star formation efficiency within the halos de- 
creases in time (for instance by feedback mechanisms) , could 
reduce the gradient dSTt/dv. By arbitrarily adjusting the ef- 
ficiency with time it is possible to perfectly mimic the slower 
DM heating, however it is unclear how to physically moti- 
vate such a rapidly declining star formation efficiency with 
time. We will explore this aspect in future work. 
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Figure 7. The concentration function adopted to calculate the 
boost-factor in Eq. 3. 



4.5 Observational perspectives 

We find that the most promising of our annihilating DM 
candidates produce a strong deviation in the HI 21 cm line 
mainly in two redshift ranges: (i) z ~ 60 — 200 (correspond- 
ing to z/ ~ 10 - 25 MHz) with a DM signal ATt.oM ~ 
10-20 mK; (ii) z < 18-30 {u ~ 45-80 MHz), where the ab- 
sorption feature produced by galaxies is strongly suppressed, 
producing a signal ATb^DM 20 — 110 mK depending on 
the DM particle. 

The first high redshift signature is at ~ 10 — 25 MHz, 
a frequency too low to be observed by ground based radio 
interferometers due to the severe distortions produced by 
terrestrial ionosphere, which at frequencies ^ 30 MHz is 
known to make any observation virtually impossible. Re- 
cently a new generation of space or Moon based radio in- 
terferometers able to probe the Dark Ages at very low fre- 
quencies has been proposed, such as the Dark Ages Lunar 
Interferometer (DALI, Lazio et al. 2007), the Lunar Radio 
Array (LRA, Lazio et al. 2009) and the Dark Ages Radio 
Explore (DARE, Burns et al. 2012). These instruments, if 
approved and built, would be located on the dark side of 
the Moon, overcoming two major hurdles of low frequency 
radio observations, i.e. radio-frequency interference (RFI) 
from human-generated radio signals, and ionospheric dis- 
tortions. They could therefore probe the Dark Ages up to 
z ~ 100. 

The second signature at 45 — 80 MHz could be 

observed by current or planned radio observatories. We focus 
here on some of the main experiments. 

The LOw Frequency ARray (LOFAR), based in the 
Netherlands, has two observational bands, a low band (30- 
85 MHz) and a high band (115-230 MHz). The high band 
array is the one that will be devoted to the EoR experiment, 
measuring the redshifted HI 21 cm radiation at 6 ^ z Si 11.4 
with a resolution of ~ 3 arcminutes on a field of view of ~ 120 
square degrees (Zaroubi 2012). The low band array probes 
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Figure 8. (Top panel) Brightness temperature STt as a func- 
tion of frequency u assuming the "fiducial" reionization model for 
three cases: (i) a standard one without DM annihilations (black 
solid line); (ii) a case including 10 GeV DM annihilations with 
a standard {a-v)tii and galactic heating and hya coupling (red 
solid line); (iii) a final scenario with 10 GeV DM annihilations 
with {cFv)ii^ assuming that galaxies produce Ly« coupling with- 
out heating the gas (red dotted line). (Bottom) Gradient dSTi^/dv 
of the same quantities. The different slopes relative to heating 
from DM or galaxies is evident 

the frequencies at which the DM signal is stronger, however 
the detection of HI 21 cm signals at such low frequencies will 
be beyond reach because of a substantial drop in sensitivity 
(see e.g. Zaroubi 2010 and references therein). 

The Murchison Widefield Array (MWA), located in the 
radio-quiet Western Australia Outback, will observe at fre- 
quencies from 80 to 300 MHz and therefore misses the fre- 
quency bands in which a DM signal could be detected. 

The Square Kilometer Array (SKA) which is planned to 
be completed in 2024 in Australia, New Zealand and South 
Africa, will probe a frequency range from 70 MHz to 10 
GHz. The antenna elements that will make up the "SKA- 
low" array will cover the frequency range = 70 — 200 MHz 
and will achieve a sensitivity an order of magnitude higher 
than previous experiments. The design is still preliminary 
and the telescope is at least a decade away from being fully 
operational, however SKA could detect a DM signal in the 
HI 21 cm line, if frequencies down to 70 MHz or lower will 
be probed. 

The main scientific goal of the aforementioned inter- 
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Table 2. DM signal for the considered DM models. 
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ferometers is however to measure the spatially fluctuating 
components of the HI 21 cm signal. On the other hand the 
detection of the angularly averaged all-sky signal could be 
performed by a single dipole, the main challenge being a 
precise enough calibration to extract the signal from the 
backgrounds. This simpler and lower cost alternative is ex- 
plored by instruments such ad EDGES and the Cosmolog- 
ical Reionization Experiment (CoRE, Chippendale 2009), 
while more advanced second generation instruments such 
as the Large-aperture Experiment to Detect the Dark Ages 
(LEDA) are under construction. A Moon orbiting space ob- 
servatory such as the Dark Ages Radio Experiment (DARE) 
seems to be a viable option in the near future, in contrast 
with the much more complex and expensive task of placing 
a large interferometer on the dark side of the Moon (DALI, 
LRA). 

The Experiment to Detect the Reionization Step 
(EDGES) located at the Murchison Radio-astronomy Ob- 
servatory in Western Australia, measured the radio spec- 
trum between 100 and 200 MHz with all systematic trends 
in the measurement reduced to below 75 mK, and allowed 
to exclude a rapid reionization timescale of Az < 0.06 at the 
95% confidence level. The EDGES team will attempt in the 
next few years to reduce the systematics of over one order of 
magnitude and to push observations up to z ~ 20 (see e.g. 
Bowman & Rogers 2010; McQuinn 2010; Pritchard & Loeb 
2010). 

One of the most promising instruments for the detection 
of a DM signal from the Dark Ages is the Large-Aperture 
Experiment to Detect the Dark Ages (LEDA), a proposed 
array that will cover a band ~ 45 — 90 MHz and that 
will feature array-based calibration to improve the accu- 
racy of foreground subtraction from the total-power signal 
(Greenhill & Bcrnardi 2012). The fact that LEDA is opti- 
mized for the detection of the all-sky HI 21 cm signal to- 
gether with its low frequency capabilities make it the ideal 
ground based instrument to detect the HI 21 cm deviations 
induced by DM annihilations. 



5 SUMMARY AND CONCLUSIONS 

We consider three of the most popular annihilating DM 
candidates as sources of X-ray, ionizing and Lya radiation: 
thanks to our detailed Monte Carlo treatment MEDEA2 we 



computed their energy deposition into the IGM as a func- 
tion of redshift and ionized fraction. We were then able to 
calculate precisely and fully consistently the effects of the 
energy release from DM annihilations on the IGM thermal 
history and, ultimately, the imprint on the HI 21 cm line 
at 20 2 ^ 200. For each DM candidate we compute with 
a modified version of the publicly available code CAMB the 
deviations produced by each DM model on the CMB power 
spectrum. This allows us to select two values for the ther- 
mally averaged annihilation cross section, the standard ther- 
mal case {av)th = 3 x 10~^^cm^ s~^ and the maximum value 
{av)max that produces deviations in the CMB power spec- 
trum within 3(T of the WMAP 7 results. We combine each 
DM model with two prescriptions for radiation from astro- 
physical sources (computed with the public code 21cmFAST 
). We investigate how the formation of luminous sources af- 
fects the imprint of DM annihilations on the IGM, and study 
ways to disentangle the DM signal in this realistic scenario. 

Finding a strong deviation in the absorption feature in 
the all-sky averaged STb at 30 ^ z ;S 200 will be a clear sign 
of energy release by DM during a cosmic phase in which no 
collapsed sources of radiation have yet formed. However it 
is even more interesting to study the joint effects of energy 
injection into the IGM from galaxies and DM annihilations 
at 20 z ^ 30, as we find that the Lya flux from the first 
luminous sources combined with the early heating from DM 
annihilations produces strong deviations in the HI 21 cm sig- 
nal at a redshift range which could be realistically probed 
by future radio observations in the next few years. Next gen- 
eration radio facilities such as DALI, LRA or DARE would 
certainly have ideal characteristics for the detection of a DM 
signal, but even current or already funded instruments such 
as MWA, LOFAR and in particular SKA and LEDA would 
achieve the sensitivity to probe a ATi,^bm ~ 20 — 110 mK 
at ~ 45 — 80 MHz, if these frequencies will be observed 
by their flnal designs and assuming that the strategies for 
foreground removal and ionospheric corrections will be suc- 
cessful. 

We summarize our results as follows. 

• DM annihilation signal. 

(i) Our least massive DM candidate, the 10 GeV Bino like 
neutralino, produces a large signal in the "fiducial" model, 
both when assuming the standard thermal cross section 
{av)th and when taking into account the maximum cross 
section allowed by CMB data, {av)max - In the first case the 
DM 21 cm signal ATb,DM is of the order of 5 — 10 mK on a 
redshift range 45 < z < 300 while ATt^DM ~ 10-100 mK at 
20 ^ 2 ;S 30, with the peak signal occurring at a frequency 

~ 80 MHz. Assuming the higher annihilation cross section 
{o'v)max increases ATt^oM and totally erases the second ab- 
sorption feature at 20 z ;$ 30, the lack of which represents 
a very strong - and detectable - DM signature. For the "ex- 
treme" reionization model we find in general a smaller sig- 
nal with a peak AT^^om 25 (50) mK at z ~ 21 assuming 

{av)th ({(Tv)max)- 

(ii) The annihilating 200 GeV Wino produces very different 
results depending on the (av) assumption of a thermal or 
maximal value. In the first case ATt,DM is negligible except 
for a weak signal ~ 15 mK a.t z ^ 23. In the second case 
instead ATi^^bm ~ 140 mK at z --^ 23 and the second ab- 
sorption feature driven by Lya coupling from stars is almost 
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completely erased. For the extreme case we find a peak sig- 
nal ATb,Di\i ~ 5(45) mK at 2 ~ 21 for {(Tv)fh {{(Tv)max)- 
(iii) The results for our most massive candidate, the heavy 
DM particle of rest mass 1 TeV that pair annihilates into 
leptons, are the least promising, with a negligible signal un- 
der the assumption of {av)th, and a ATi,^r>M ~ 60(20) mK 
at z ~ 23(100). For the extreme case we find a peak signal 
ATb^DM ^ 0.5 (10) mK at z ~ 21 for (av)th {{(^v)max)- 
We summarize our results in Table 2. 

• Strategies to isolate the DM signal. 
We find that some of our results for a specific DM model 
working in the framework of the "fiducial" reionization 
model (e.g. the case of 10 GeV annihilations assuming 
{o'v)th) are hard to distinguish from the standard signal 
STbfi assuming the "extreme" reionization history. However 
when we study the heating rates per baryon for 10 GeV 
DM annihilations (edm) and for galactic X-ray heating (ex) 
we notice that the heating rate from astrophysical sources 
of X-rays evolves much quicker than the DM heating, since 
the fractional increase of the collapsed fraction in ^ IMq 
halos which drive the DM heating is much slower than the 
fractional increase of the halos which host the first galax- 
ies. This crucial difference presents us with a clean way 
of disentangling heating by DM annihilations from heat- 
ing by galaxies. To better quantify this point we study the 
case of annihilating 10 GeV Binos with {av)th, keeping the 
Lya pumping due to stellar sources but switching-off astro- 
physical X-ray heating, which could indeed take place at 
lower redshifts without breaking any current observational 
constraints. By coincidence STi, starts increasing at z ^ 22 
{u ~ 60 MHz), same as the case in which source heating 
is taken into account, however DM heating happens with 
a different slope. We therefore study the gradient of the 
signal as a function of frequency, dSTb/dv, and find that 
when neglecting DM heating the slope reaches values of 
dSTb/dv ~ 14 mK/MHz, while in the case in which DM 
annihilations are taken into account the slope is confined to 
a considerably smaller dSTb/dv ^ 4 mK/MHz. Although it 
is possible in principle to mimic the slower DM heating by 
arbitrarily reducing the star formation efficiency with time, 
it is unclear whether such an "ad-hoc" model could be physi- 
cally motivated, and we defer this for future analysis. We can 
therefore conclude that the detection of a small peak gra- 
dient (dSTb/dv ~ few mK/MHz) will be a strong evidence 
of DM annihilation, an exciting prospect that could become 
reality in the next few years thanks to a new generation of 
radio instruments such as LEDA and DARE, optimized for 
the detection of the global HI 21 cm signal during the Dark 
Ages. 
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